J Wood Sci (2007) 53:427-435 
DOI 10.1007/sl0086-007-0888-6 


© The Japan Wood Research Society 2007 


ORIGINAL ARTICLE 


Shinya Sasaki • Motoyuki Shimizu • Hiroyuki Wariishi 
Yuji Tsutsumi • Ryuichiro Kondo 

Transcriptional and translational analyses of poplar anionic peroxidase 
isoenzymes 


Received: November 15, 2006 / Accepted: January 22, 2007 / Published online: May 29, 2007 


Abstract Anionic peroxidases have been proposed to be a 
key enzyme for lignification in poplar. On the other hand, 
there are many genes encoding an anionic peroxidase in 
Populus trichocarpa genome, and their physiological func¬ 
tions are still unknown. Ampholine isoelectric focusing 
analysis revealed anionic peroxidases as dominant peroxi¬ 
dases in enzyme preparations from various organs. Using 
two-dimensional electrophoresis (2-DE) followed by pep¬ 
tide mass fingerprint (PMF) analysis, we surveyed the local¬ 
ization of anionic peroxidase isoenzymes in various organs 
of Populus alba L. Peroxidase isoenzymes were extracted 
from various organs and fractionated by a Concanavallin A 
Sepharose column. Each protein was separated by 2-DE 
gels and some anionic peroxidase isoenzymes in each organ 
were identified via PMF analysis. Transcript and protein of 
individual peroxidase indicate that the expression profile of 
each isoenzyme is quite different, for example, organ- 
specific gene, stress-response gene, and multifunction gene, 
even though they are in the same cluster. These results sug¬ 
gest that individual anionic isoenzymes in this small cluster 
were differently regulated at transcription, translation, or 
posttranslation. 

Key words Anionic peroxidase • Lignification • Poplar • 
Real-time PCR • Two-dimensional gel electrophoresis 


Introduction 

One of the most studied functions of plant peroxidases is 
their role in lignification. Extensive studies, so far, have 
been carried out to identify whether some specific peroxi- 
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dase isoenzymes are correlated with lignification in many 
plants, especially tobacco and poplar. “ In these studies, 
anionic peroxidases have been discussed as being responsi¬ 
ble for lignification because they have been often found as 
highly transcribed peroxidase genes in the developing xy- 
lem or as dominant peroxidase proteins in the enzyme 
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preparation from the same organ. “ One of anionic peroxi¬ 
dases, Populus kitakamiensis prxA3a gene, has been re¬ 
ported as a highly transcribed and translated gene in the 
developing xylem. 3,4 Li et al. 3 generated transgenic poplar 
trees {Populus sieboldii x Populus gradidentata ) by an 
antisense strategy with a prxA3a gene and achieved a 20% 
reduction in lignin content of the transgenic plant at maxi¬ 
mum when compared with control trees. These studies have 
proposed that anionic peroxidases are key enzymes in 
lignification. 36 

On the other hand, it has been known for a long time 
that plant peroxidases exist as a large family of isoenzymes 
in plants; for instance, 73 peroxidases and 138 peroxidase 
genes have been identified throughout the Arabidopsis ge¬ 
nome and the Oryza sativa genome, respectively. 7,8 Very 
recently, the genomic sequence of Populus trichocarpa was 
completed by the United States Department of Energy 
(USDOE) Joint Genome Institute, 9 and the preliminary 
annotation data are available to the public (http://genome. 
jgi-psf.org/Poptrl_l/Poptrl_l.home.html). In our laborato¬ 
ry, at least 82 sequences encoding complete peroxidase pro¬ 
tein were found with manual curation within the database 
(unpublished data). It was also found that 8 genes were 
highly similar to prxA3a among them. Therefore, it is of 
interest to investigate whether physiological functions of 
these isoenzymes are the same as prxA3a , because protein 
structure and substrate specificity of these isoenzymes are 
presumed to be highly similar. 

Recently, some studies on Arabidopsis thaliana and Ory¬ 
za sativa peroxidases were performed using the microarray 
and reverse-transcribed polymerase chain reaction (PCR) 

n 4 A 1 1 

technique. ’ ’ High density microarray analyses provide 
qualitative and semiquantitative data for large numbers of 
mRNAs; however, the levels of gene expression do not 
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always directly reflect active enzyme levels. Enzyme activ- 
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ity is often regulated through posttranslational modifica¬ 
tions such as phosphorylation and acylation or through 
protein sorting, protein-protein interactions, and/or con¬ 
trolled proteolysis. 14 This is also true of plant peroxidases. 10 
Therefore, the double expression analyses of proteins and 
transcripts will afford better understanding of the function 
of peroxidase gene expression. 

In this study, we investigated protein and transcript 
level of the individual anionic peroxidase isoenzyme that 
is highly similar to prxA3a in various organs in the poplar 
tree to estimate the physiological function of these 
peroxidases. 


Materials and methods 

Plant materials 

Poplar ( Populus alba L.) callus was induced and developed 
on Murashige and Skoog basal medium supplemented 
with 3% sucrose, 1.0 ppm 2,4-dichlorophenoxyacetic acid, 
0.5 ppm kinetin, and 0.8% agar. The callus was maintained 
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on the medium at 25°C in darkness. The poplar tree 
(P. alba) grown at Kyushu University was used for experi¬ 
ments of peroxidase gene cloning, transcription, and two- 
dimensional electrophoresis (2-DE) analyses. 

Cloning of poplar peroxidases 

Total RNAs from xylem, shoot, leaf, bark, petiole, and 
stress-treated leaves were prepared by Concert Plant RNA 
Reagent (Invitrogen, Carlsbad, CA, USA) according to the 
manufacturer’s instructions. cDNA was synthesized in a fi¬ 
nal volume of 20 pi that included 1 pg of total RNA, 1 pM 
oligo-dT 18-mer primer, 10 units of RNase inhibitor, and 10 
units of AMV Reverse Transcriptase (Takara, Japan) ac¬ 
cording to the manufacturer’s instructions. The PCR prim¬ 
ers were selected in the sequence region specific to each 
peroxidase homolog of Populus trichocarpa. The following 
oligonucleotide primers were used in the PCR reaction: 
PAPOl and PAP02 sense primer (5'-GAAAAATG 
CATCTTTCTA AGGC-3'), PAPOl and PAP02 antisense 
primer (5'-TGTC A AC A AG ACTCC ATTTGG-3'), PAP 03 
sense primer (5'-CTGGTTTGGAGGATCACTTCC-3'), 
PAP03 antisense primer (5'-TTT GAT CT C A AT CTT GT C 
TTCTC-3'), PAP04 sense primer (5'-GAGAAAATG 
AAGCTTTCAAAAC-3'), PAP04 antisense primer (5'- 
C A ATTTTTTCTTGGC ATAGTCTGGC-3'), PAP 05 
sense primer (5'-GTTCTGTTCA AG A AG A A A A ACA AT 
GTC-3'), PAP05 antisense primer (5'-CTATCTAGTG 
ATCGCATAGC-3'), PAP06 sense primer (5'-CTCATAG 
CATAGAATTTCCGAAGACC-3'), and PAP06 anti- 
sense primer (5'-AT GT CTT CT GTT CT AGCT AC AGT-3'). 
All PCR amplifications were carried out in 30 pi of reaction 
mixture containing primers (0.2 pM), dNTPs (0.2 mM each), 
and 1 unit of Takara Ex-Taq polymerase. 

The PCR procedure started with 5 min of denaturation 
at 94°C and was carried out over 35 cycles of 30 s of dena¬ 


turation at 94°C, 30 s of annealing at 55°C, and 1.5 min of 
extension at 72°C. The PCR products were purified from 
the agarose gel with a GEL-M kit (Viogen, Japan) accord¬ 
ing to the manufacturer’s instructions and were subcloned 
into pXcmknl2 vector. The resulting ligation product was 
transformed into Escherichia coli strain DH5a according to 
the manufacturer’s protocols. The clones were sequenced 
by a dideoxy dye termination method (Thermo Sequence 
Cycle Sequence Kit, Amersham Bioscience, Piscataway, 
NJ, USA) with a LIC-4000 sequencer (Aloka, Japan). The 
region representing the signal peptide of the peroxidase 
homolog was predicted with a primary structure analysis 
program, SignalP (http://www.cbs.dtu.dk/services/SignalP/). 
The alignment was created with the GENETYX-Mac 
version 13. 0. 3. (Genetyx, Japan). The name of each se¬ 
quence corresponds to the deduced protein of the following 
DNA accession numbers: PAPOl: AB206042, PAP02: 
AB260988, PAP03: AB260989, PAP04: AB260990. 
PAP05: AB260991, PAP06: AB260992. 

Stress treatment 

Leaflets were cut off and subjected to wounding, cadmium, 
and H 2 0 2 treatments. For wounding treatment, fully ex¬ 
panded leaves of healthy poplar plants were detached and 
immediately wounded with a needle-point holder. Wound¬ 
ed leaves were put on wet paper towel in petri dishes and 
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incubated at 25°C in darkness. For cadmium treatment, 
leaves were submerged in 11 of aqueous 1 mM CdCl 2 solu¬ 
tion at 25°C in darkness. For water treatment, leaves were 
submerged in 11 of water at 25°C in darkness. For H 2 0 2 
treatment, leaves were submerged in 11 of aqueous 200 mM 
H 2 0 2 at 25°C in darkness. All stress treatments were con¬ 
ducted for 1, 5, and 10 h, respectively, and the leaflets were 
then frozen in liquid nitrogen. For the control, leaflets, in¬ 
ner bark, petioles, shoots, and xylem were harvested, then 
immediately frozen in liquid nitrogen. In the preparation of 
xylem, bark was removed from the stem and thin flakes 
(thickness <lmm) were whittled from the surface of the 
xylem using a feather knife. 

Preparation of glycoproteins 

Xylem, petiole, leaf, and inner bark organ were ground with 
liquid nitrogen and homogenized in 50 mM Tris-HCl buffer 
(pH 7.5). After centrifugation at 20000g at 4°C for 30min, 
the supernatant was salted out with ammonium sulfate with 
80% (w/v) saturation. The protein precipitate was collected 
by centrifugation (20 000 g for 30 min), then dissolved in 
50mM Tris-HCl buffer (pH 7.5). After centrifugation at 
40 000g at 4°C for 30 min, the supernatant was desalted with 
a PD-10 column (Amersham Bioscience). The glycopro¬ 
teins were separated by a Concanavallin A Sepharose col¬ 
umn with the 50 mM Tris-HCl buffer (pH 7.5) containing 
0.5 M NaCl and 0.2 M tz-methylglycopyranoside. Finally, 
purified glycoproteins were desalted by a PD-10 column 
with 50 mM Tris-HCl buffer (pH 7.5) and concentrated via 
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ultrafiltration (molecular weight cutoff <10000, Advantec, 
Japan). 

The solution was heated at 80°C for 5 min and insoluble 
material was removed by centrifugation (15 000g for 30 min). 
Four volumes of cold acetone (-20°C) were added and the 
solution was incubated overnight (-20°C). After centrifuga¬ 
tion (15000g for 30min), the precipitate was washed with 
cold acetone (-20°C), and the pellet was solubilized in 2- 
DE solubilization buffer containing 7M urea, 2M thiourea, 
4% CHAPS, 2% DTT, 0.5% IPG buffer (pH 4.0-7.0, pH 
3.0-5.0; Amersham Bioscience), and a trace of bromophe- 
nol blue. The sample was then incubated for 2h at room 
temperature, and insoluble material was removed by cen¬ 
trifugation (15000g for 30min). 


Ampholine isoelectric focusing analysis 

Isoelectric focusing (IEF) gel (pH 3.5-9.5) electrophoresis 
was performed according to the manufacturer’s instruc¬ 
tions. The gel was mounted on a water-cooled (10°C) elec¬ 
trophoresis unit (2117 Multiphor II; Amersham Bioscience). 
Protein samples were subjected to electrophoresis at 1500 V, 
30 W, and 50 mA for 90 min. Then the gel was stained with 
10 mM guaiacol and ImM hydrogen peroxide or coomassie 
brilliant blue R-250 reagent at room temperature. Peroxi¬ 
dases were visualized and the isoelectric points were esti¬ 
mated by reference to the mobility of marker proteins 
(Amersham Bioscience). Peroxidase activity was assayed 
using guaiacol as a substrate according to described 
methods. 18 


Two-dimensional gel electrophoresis 

Two-dimensional gel electrophoresis was performed as pre¬ 
viously described. 19 Isoelectric focusing was carried out with 
an IPGphor system (Amersham Bioscience). Immobilized 
pH gradient strips (pH 4.0-7.0, pH 3.0-5.6, pH 3.5-4.5, 
18 cm; Amersham Bioscience) were rehydrated for 12 h, and 
then 200 pg of protein in 2-DE solubilization buffer was fo¬ 
cused in four steps at 500 V (lh), 500-1000V (lh), 1000- 
8000 V (2h), and 8000 V (8h). After completion of focusing, 
strips were equilibrated with buffer containing 6M urea, 
130 mM DTT, 30% glycerol, 2% sodium dodecyl sulfate 
(SDS), and a trace of bromophenol blue then with a buffer 
containing 6M urea, 135mM iodoacetamide, 30% glycerol, 
2% SDS, and a trace of bromophenol blue. The strips were 
loaded onto precast 10% or 12% homogenous polyacryl¬ 
amide gels (20 x 20 cm). The lower running buffer contained 
385mM Tris, 50mM glycine, 0.1% SDS, and 0.02% sodium 
azide, while the upper running buffer was identical except 
it lacked sodium azide. The system was run at 1000 V and 
24 mA per gel. Gel slabs were stained in 7.5% acetic acid 
solution with 0.02% SYPRO Red (Takara) and incubated 
with gentle rocking at room temperature for 1 h. After re¬ 
moval of the staining solution, gels were washed in 7.5% 
acetic acid solution for 30 min. 


In-gel tryptic digestion 

In-gel tryptic digestion was performed as previously de¬ 
scribed 20 with a slight modification. The target spot was ex¬ 
cised and cut into 2-mm cubes and the gel pieces were 
transferred into a 1.5-ml microcentrifuge tube and washed 
with 40% 1-propanol at room temperature for 15 min. After 
removal of 1-propanol solvent, 200 mM ammonium bicar¬ 
bonate in 50% acetonitrile was added, and the sample was 
incubated at room temperature for 15 min. The gel pieces 
were then dried and covered with 20 ng/pl modified trypsin 
(Promega, Madison, WI, USA) in a minimal volume of 
100 mM ammonium bicarbonate to rehydrate the pieces; 
they were then incubated at 37°C. After a 12-h incuba¬ 
tion, the supernatant was collected, and the gel pieces 
were extracted once with 100 mM ammonium bicarbo¬ 
nate, followed by two extractions with 80% acetonitrile 
containing 0.05% trifluoroacetic acid. The supernatant and 
extracts were combined and concentrated to the required 
concentrations. 

Matrix assisted laser desorption ionization time-of-flight 
mass spectrometry (MALDI-TOF-MS) analysis 

The resulting peptide mixtures were desalted using ZipTips 
C18 (Millipore, Billerica, MA, USA) and eluted onto a 96- 
well MALDI target plate. Samples of 2 pi were mixed on 
the plate with lpl of 10 mg/ml tz-cyano-4-hydroxycinnamic 
acid solution in 0.1% trifluoroacetic acid in H 2 0/acetoni- 
trile (1:1). Samples were then dried at room temperature. 
Mass spectral data were obtained using a Voyager DE mass 
spectrometer equipped with a 337 nm N 2 laser in the posi¬ 
tive ion reflectron mode (Applied Biosystems, Foster City, 
CA, USA). Spectral data were obtained by averaging 64 
spectra, each of which was the composite of 64 laser firings. 
Internal mass calibration was performed using bradykinin 
(904.45 Da) and adrenocorticotropic hormone (ACTH) 
(2465.75 Da). 

Construction of putative in silico peroxidase 
protein library 

The preliminary genomic annotations of Populus tricho- 
carpa were opened to the scientific community by the 
USDOE Joint Genome Institute (JGI) (http://genome.jgi- 
psf.org/Poptrl_l/Poptrl_l.home.html). 9 One hundred and 
eight peroxidase genes were listed as candidates by the JGI 
site. The lists of amino acid sequences of peroxidase genes 
and other genes were extracted as a fasta file, after which 
the data were converted into a MASCOT file. 

Identification of proteins 

Peptide mass fingerprint (PMF) was utilized for protein 
identification by analyzing the size of tryptic fragments via 
the MASCOT search engines (Matrix Science, Boston, MA, 
USA) using the entire NCBI protein database and the Pop- 
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ulus trichocarpa in silico protein library. For effective PMF 
analysis, it was assumed that peptides were monoisotopic, 
and the possibility that methionine residues were oxidized 
was considered. The fingerprinting method allowed for a 
maximum of one missed tryptic cleavage per protein. The 
maximum deviation permitted for matching the peptide 
mass values was 100 ppm. Scores greater than 65 were con¬ 
sidered to be significant. 


Real-time PCR analysis 

Total RNA extraction and reverse-transcriptase reaction 
were performed as described in the cloning of poplar per¬ 
oxidases. Real-time PCR was performed in a final volume 
of 10 pi with a Line Gene (Bio Flux, Japan). The SYBR 
Premix Ex Taq (Takara) was used according to the manu¬ 
facturer’s instructions with a final concentration of 0.2 pM 
for each primer. PCR amplification was performed as fol¬ 
lows: (1) an initial denaturation at 95°C for 5 min, (2) 45 
cycles, with 1 cycle consisting of denaturation 95°C for 15 s, 
annealing at 55°C for 15 s, and elongation at 72°C for 30 s. 
Amplicon specificity was verified by melting-curve analysis 
conducted at 65° to 95°C with stepwise fluorescence acquisi¬ 
tion and by 2% agarose gel electrophoresis staining with 
ethidium bromide. No fluorescence was detected from real¬ 
time PCR amplification without a template. The PCR prim¬ 
ers were selected in the sequence region specific to each 
peroxidase homolog and the PCR product was designed at 
a length of 240-258 bp. The following oligonucleotide prim¬ 
ers were used for real-time PCR analysis: each primer of 
PAPOl and PAP02 was 5'- GTT GTT GAT AG A AT G A A 
GGCTTTG-3' (forward) and 5'-GCTGACATTAGTGAA 
GCTCTC-3' (reverse), respectively; each primer of PAP03 
was 5'-GATAACATGAAGGCTGCAG (forward) and 
5'-TACTGGTGTTGAGGCCCACGGCT-3' (reverse), re¬ 
spectively; each primer of PAP04 was 5'-GTTGATATA 
ATGAAGGAAAGG-3' (forward) and 5'-GTTGATATA 
ATGAAGGAAAGG-3' (reverse), respectively; each 
primer of PAP05 was 5'-GTCGAACATTTAGCAACA 
GGC-3' (forward) and 5'-T G ACC A A AGT G ACCGTGG 
CAGCG-3' (reverse), respectively; each primer of PAP06 
was 5'-CAAACCCCGATCCAACCTTAA-3 / (forward) 
and 5'-CTATCGAAGAAGGCGGCCTGG-3 / (reverse), 
respectively. ACTIN was used as a reference gene (acces¬ 
sion number AB025795). Each primer sequence was 5'-GC 
CCAGAGGTCCTCTTCCAA-3' (forward) and 5'-GGGG 
CTAGTGCTGAGATTTCCTTGC-3' (reverse), respec¬ 
tively. Each reaction was repeated three times. The data 
were normalized by actin transcripts that were assumed to 
be constantly expressed in all organs, and the relative ex¬ 
pression of the transcript in various stress treatments was 
calculated. Template DNA for a working curve was synthe¬ 
sized by the PCR technique. All PCR amplifications were 
carried out in 30 pi of reaction mixture consisting of each of 
the primers (0.2pM), dNTPs (0.2mM each), and 1 unit of 
Takara Ex-Taq polymerase. The synthesized PCR product 
was electrophoresed by 2% agarose gel staining with ethid¬ 
ium bromide. The PCR product was purified from the aga¬ 


rose gel with a GEL-M kit according to the manufacturer’s 
instructions. The concentration of the PCR product was 
measured by spectrophotometer and used to establish a 
working curve. 


Results and discussion 

Peroxidase isoenzymes in enzyme preparations from 
various organs 

When 2-DE followed by PMF analysis of the crude extract 
from poplar developing xylem was performed, over 20 pro¬ 
teins, such as enolase, triose-phosphate isomerase, and al¬ 
cohol dehydrogenase were tentatively identified. However, 
no peroxidase protein was found among them (data not 
shown). In the Ampholine IEF followed by active staining, 
anionic peroxidases were found as dominant peroxidase 
proteins in the enzyme preparation from various organs 
(Fig. 1). In the inner bark organ, only anionic peroxidases 
were observed. Indeed, anionic peroxidases (pi 4.0-5.0) 
were also dominantly observed in other organs (Fig. 1). 
These results indicated that expression of peroxidase pro¬ 
teins was lower than that of major proteins (for instance, 
glycolytic enzymes and primary metabolism), and it appears 
that partial purification of peroxidase isoenzymes from 
crude protein preparation is needed in order to improve the 
detection limit for 2-DE analysis. It is known that plant 
peroxidases are glycoproteins. Therefore, peroxidase isoen¬ 
zymes extracted from leaf, bark, xylem, and petiole were 
fractionated by a Concanavallin A Sepharose column. This 


Fig. 1 . Isoelectric focusing 
(IEF) analysis of crude 
proteins from various organs 
of Populus alba. Crude 
proteins extracted from 
various organs (guaiacol 
oxidation activity: 0.06 unit) 
were subjected to IEF 
analysis. After electrophore¬ 
sis, the gel was stained with 
10 mM guaiacol and ImM 
hydrogen peroxide at room 
temperature. Individual 
bands correspond to 
peroxidase isoenzymes 
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Table 1. Relative peroxidase activity of crude enzyme and Con A purified fraction (%) 


Fraction 

Leaf 

Petiole 

Bark 

Xylem 

Crude fraction 

100 

100 

100 

100 

ConA purified fraction 

90 ±2 

85 ±5 

91 ±8 

63 ±5 

ConA flow through fraction 

ND 

ND 

ND 

30 ±3 

Yield of protein 

12.4 

11.5 

9.8 

15.4 


The peroxidase activity of each crude enzyme is defined as 100% 
ND, not detected 


Table 2. Peroxidases identified in various organs 


Spot no. 

Protein name 

Accession number 
{P. alba) 

JGI annotated gene (P. trichocarpa) 

tpl 

tMw 

Coverage 3 (%) 

Score 

1 

PAPOl 

AB206042 


4.03 

33 944.36 

38 

71 

2 

PAP02 

AB260988 


4.06 

33 885.35 

35 

68 

3 

PAP03 

AB260989 


4.51 

33 572.05 

23 

90 

4 

Calreticulin 


estExt_Genewisel_vl.C_LG_XIII0633 

4.35 

47 425.07 

29 

103 


a Sequence coverage in peptide mass fingerprint analysis 


separation step gave a concentration effect of 7.3 to 9.4 
times and 90% recovery of peroxidase activity except that 
ca. 30% of activity in xylem did not adsorb onto the column 
(Table 1). However, it was confirmed that the main isoen¬ 
zymes from xylem were recovered in the purified fraction 
in the Ampholine IEF analysis (data not shown). The frac¬ 
tionation by using a Concanavallin A Sepharose column 
appears to be appropriate for the concentration of peroxi¬ 
dases from crude enzyme preparations as a prefractionation 
procedure. 

In Fig. 2a, high expression of spots nos. 1 and 2 were 
detected in callus. Individual spots were subjected to 
PMF analysis. These spots were most similar to EU- 
GENE3.00010103 or GRAIL3.0116010201 (http://genome. 
jgi-psf.org/Poptrl_l/Poptrl_l.home.html) of Populus 
trichocarpa peroxidases. However, the score of PMF was still 
low via a MASCOT search. Therefore, cDNAs of these 
peroxidase isoenzymes of Populus alba were cloned. 
The cluster containing EUGENE3.00010103 and GRAIL 
3.0116010201 consisted of anionic peroxidase genes 
(EUGENE3.00010103, GRAIL3.0116010201, ESTEXT_ 
GENE WISE1_V 1. C_LG_III2714, ESTEXT_GENE- 

WISE1_V1.C_LG_III2066, ESTEXT_GENEWISE1_V1. 
C_LG_XVI1838, EUGENE3.00031890, ESTEXT_ 
F GENESH1_PG_V 1. C_LG_I0098, and EUGENE 
3.00010106). Based on the their DNA sequence we designed 
specific primers for each peroxidase homolog of P. tricho¬ 
carpa , and six cDNAs encoding mature anionic peroxidases 
were cloned from P. alba. 

Again, the PMF data of the protein spots were subjected 
to the MASCOT search of our peroxidase protein library 
including the cloned P. alba peroxidases. Spots nos. 1 and 
2 were identified as the PAPOl protein and PAP02 pro¬ 
tein, respectively, by PMF analysis (Table 2). Although 
PAPOl and PAP02 proteins were highly similar to each 
other (97.3% homology at amino acid level), MALDI- 
TOF-MS analysis was able to clearly detect the different 
molecular weights of the tryptic peptides (PAPOl: 
150-DAANAFLEPPTLTLDQLR-167, M w = 1969.25; and 


PAP02: 150-DAANAFL/ , PPTLTLDQLR-167, M w = 
1953.20) from two peroxidases, respectively. In bark and 
petiole, the 2-DE profiles were nearly equivalent (Fig. 2b, 
c), and PAPOl and PAP02 proteins were identified 
(Table 2). In the leaf, the 2-DE profile was slightly different 
from the profiles of bark and petiole (Fig. 2b, c, d), and the 
PAPOl protein was detected but not the PAP02 protein. 
The 2-DE profile of xylem (Fig. 2e) was very different from 
those of other organs. The PAPOl protein was also de¬ 
tected in the xylem, but its expression level was obviously 
lower than in the other organs (Fig. 2b, c, e). The PAPOl 
protein was expressed in all organs. On the other hand, the 
PAP02 protein was observed in the bark and petiole. 
PAPOl and PAP02 were highly similar to the PXP3-4 in 
P. trichocarpa [accession number X97350, 94.6% homology 
at nucleotide level {PAPOl vs PXP3-4 ), 94.3% homology 
at nucleotide level {PAP02 vs PXP3-4)\. PXP3-4 was re¬ 
ported as syringaldazine-oxidizing peroxidase among pop¬ 
lar xylem peroxidases. 5,6 The transcript of PXP3-4 in xylem 
and bark is highest among other organs. 6 The high tran¬ 
scription in bark was consistent with the expression profile 
of PAP02. PAPOl and PAP02 were also highly similar to 
the PCY 2-6 gene in Populus nigra [accession number 
D83224, 97% homology at nucleotide level {PAPOl vs 
PCY2-6 ), 96.7% homology at nucleotide level {PAP02 vs 
PCY2-6 )]. Kato et al. 21 investigated the transcriptional lo¬ 
calization of PCY2-6 in P. nigra xylem and phloem by in 
situ hybridization. The mRNA was observed in the phloem 
and inner bark, and the researchers speculated that PCY2-6 
might be involved in bark development. In the 2-DE anal¬ 
ysis, PAP02 protein was observed as a highly expressed 
peroxidase in the bark and petiole compared with 
other organs. This observation is consistent with the study 
of Kato et al. Therefore, PAP02 might be involved in 
suberization. 

In the 2-DE analysis, the PAP03 protein was detected 
only in the xylem (Fig. 2e, Table 2), suggesting that PAP03 
might be one of the lignification-specific peroxidase isoen¬ 
zymes. The result of NCBI blast search showed that PAP03 
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Fig. 2. a-e. Two-dimensional 
electrophoresis of glycoproteins 
of various organs of Populus 
alba, a Glycoproteins of 
callus ( a-1 , 200 pg; u-2, 50 pg). 
b Glycoproteins of petiole 
(b-1, 200 pg; b-2, inset from b-1). 
c Glycoproteins of inner bark 
(c-1, 200 pg; c-2, inset from c-1). 
d Glycoproteins of leaf (200 pg). 
e Glycoproteins of xylem ( e-1, 
200pg; e-2, inset from e-1). 
Proteins (Nos. 1-4) were 
identified by peptide mass 
fingerprint analyses (see 
Table 2) 
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is highly similar to the prxA3a gene from Populus kinaka- 
minesis (accession number D38050, 93.9% homology at 
nucleotide level). Li et al. 3 generated transgenic poplar 
trees (Populus sieboldii x Populus gradidentata) by an an¬ 
tisense strategy with a prxA3a gene under the control of the 
original promoter. The peroxidase activities were lower in 
transgenic lines than in the control plants. In addition, lower 
lignin content and modified lignin composition was associ- 
ated with the reduction of peroxidase activity. Li et al. also 
discussed that prxA3a would be involved in the lignification 


of xylem. Takeuchi et al. 22 investigated the localization of 
the PrxA3a protein in differentiating poplar xylem, and 
found that the protein localizes in the cells in the early stage 
of secondary wall formation. Thus, the peroxidase encoded 
by prxA3a was proposed to be involved in the early stage 
of lignification. 22 NCBI blast search showed that PAP03 is 
highly similar to PXP1 in P. trichocarpa (accession number 
X97348, 92.6% homology at nucleotide level). 6 PXP1 pro¬ 
tein was found as the dominant protein in enzyme prepara¬ 
tions from xylem, 5 which is consistent with our study (Fig. 
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2e). The physiological function of PXP1 was not discussed 
in the previous reports; 5,6 however, these three highly simi¬ 
lar genes, prxA3a , PXP1 , and PAP03 , appear to be in¬ 
volved in lignification. 

In the xylem, spot no. 4 protein was identified as ca- 
lreticulin. The plant calreticulin is often ^-glycosylated. 23 
Calreticulin participates in the folding newly synthesized 
proteins and glycoproteins. It is an important component of 
the calreticulin/calnexin cycle and the quality control path¬ 
ways in the endoplasmic reticula. 

Transcript profile of identified peroxidase genes 

Three peroxidase isoenzymes were identified by 2-DE 
followed by PMF analysis and they belong to a small cluster, 
although the cluster contains at least six peroxidase genes. 
It is, therefore, of interest that the transcriptional regulation 
of the peroxidase genes in this cluster and transcription 
of these peroxidases in various organs and under the 
several stress conditions were investigated by real-time 
PCR analysis. Because PAPOl and PAP02 are very simi¬ 
lar, specific primers were not able to be designed to amplify 
these transcripts separately. Thus, consensus primers for 
PAPOl and PAP02 were designed simultaneously. Figure 
3 shows the transcription profiles of six peroxidases in this 
cluster. 

For PAPOl and PAP02 , real-time PCR analysis indi¬ 
cated that the transcript in the bark and petiole was the 
highest among various organs. The PAPOl and PAP02 
proteins were also detected as the major peroxidases in the 
bark and petiole (Fig. 2b, c). The transcript and translation 
of these peroxidase genes are highly consistent concerning 
their organ specificity. The transcript level of PAP03 was 


constant in various organs and lower than PAPOl and 
PAP02 (Fig. 3); however, the PAP03 protein was a major 
peroxidase detected only in the xylem. The transcription 
profile of PAP04 was quite similar to that of PAP03 (Fig. 
3), but PAP04 protein was not detected in xylem (Fig. 2). 
These results suggested that the expression of PAP03 and 
PAP04 appears to be controlled in the translation level or 
posttranslation level rather than in the transcript level. 

Besides lignification, many reports have suggested that 
plant peroxidases play a role in resistance to injury, patho¬ 
gens, and oxidative stress. 24 31 Therefore, the transcription 
of these peroxidase genes was determined in the stress- 
treated leaf in order to identify which peroxidase genes are 
involved in the plant defense mechanism. The transcripts of 
PAP03 and PAP04 were decreased or unchanged by all 
stress treatments (Fig. 4b, c). Under the nonstressed condi¬ 
tion, the PAP05 and PAP06 proteins were not detected 
by 2-DE analyses (Fig. 2) and their transcripts were quite 
low in real-time PCR analysis (Fig. 3). On the contrary, 
PAP05 and PAP06 were upregulated by wound stress or 
water stress (Fig. 4d, e). PAP05 was specifically upregu¬ 
lated by wounding, suggesting that this peroxidase is in¬ 
volved in the plant defense mechanism against injury. The 
transcripts of PAPOl and PAP02 were increased by 
hydrogen peroxide- and water-stress treatments (Fig. 4a). 
Figure 3 and 2-DE analysis show that PAPOl and PAP02 
were highly expressed in inner bark (Figs. 2 and 3). There¬ 
fore, PAPOl and PAP02 seem to be involved in suberiza- 
tion and the plant defense mechanism. 

The transcriptional and translational profiles of paralo- 
gous genes encoding anionic peroxidases are summarized 
in Fig. 5, which indicates that the expression of each per¬ 
oxidase in this cluster is individually regulated. Double ex¬ 
pression analyses identified organ-specific expression genes 


Fig. 3. The transcription profile 
of anionic peroxidase genes in 
various organs of Populus alba. 
The transcripts of the anionic 
peroxidase genes were quantita¬ 
tively determined by real-time 
polymerase chain reaction 
(PCR). Data are given as mean 
±SD (n = 3). The concentrations 
of individual DNA were 
calculated from a working curve 
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Fig. 4a-e. The transcription 
profiles of anionic peroxidase 
genes in stress-treated leaves of 
Populus alba, a PAPOl and 
PAP02 , b PAP03, c PAP04 , 
d PAP05 , e PAP06. The 
transcripts of the anionic 
peroxidase genes were quantita¬ 
tively determined by real-time 
PCR. The data were normalized 
by actin transcript, which was 
assumed to be constantly 
expressed in this experiment, 
and the relative expression was 
calculated for each gene in 
stress-treated leaves. Data are 
given as mean ±SD (n = 3) 
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Fig. 5. Characteristics of the expression of six anionic peroxidases in 
Populus alba. PAPOl, PAP02, and PAP03 proteins were identified 
by PMF analyses in P. alba. PAP03 and PAP04 were constitutively 
expressed, and the other genes were stress-enhanced. PAPOl and 


PAP02 genes were highly transcribed and upregulated by hydrogen 
peroxide stress and water stress. The scale bar represents 0.05 substitu¬ 
tions per amino acid residues 
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( PAP03 ), stress-response genes ( PAP05 , PAP06 ), and 
multifunctional genes ( PAPOl , PAP02) (Fig. 5). In this 
study, double expression analysis showed that positive cor¬ 
relations between transcript and their protein of peroxi¬ 
dases were not necessary. It is interesting that anionic 
peroxidases may have physiological diversity by means of 
acquiring different regulation of gene expression rather 
than modifying their protein structure and substrate 
specificity. 
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